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Abstract. We study theoretically the dynamical rectification of a terahertz AC 
electric field, i.e. the DC current and voltage response to the incident radiation, in 
strongly coupled semiconductor superlattices. We address the problem of stability 
against electric held domains: A spontaneous DC voltage is known to appear exactly 
for parameters for which a spatially homogeneous electron distribution is unstable. We 
show that by applying a weak direct current bias the rectifier can be switched from 
a zero DC voltage state to one with a hnite voltage in full absence of domains. The 
switching occurs near the conditions of dynamical symmetry breaking of an unbiased 
Q . semiconductor superlattice. Therefore our scheme allows for the generation of DC 

voltages that would otherwise be unreachable due to domain instabilities. Furthermore, 
for realistic, highly doped wide miniband superlattices at room temperature the 
generated DC field can be nearly quantized, that is, be approximately proportional 
to an integer multiple of huj/ea where a is the superlattice period and w is the AC field 
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. frequency. 
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1. Introduction 

Semiconductor superlattices (ssls) have been the subject of much research ever since 
Esaki and Tsu[l] realized that these nanostructures could exhibit B loch-oscillations [2] 
at moderate electric field strengths. This oscillatory response makes SSLs inherently 
nonlinear, and the variety of transport phenomena that has been discovered has lead 
to the point that they can now be called model systems for nonlinear transport [3] . 
Semiconductor superlattices hold great promise as amplifiers and detectors of THz 
radiation. A considerable amount of research has been put into the problem of 
THz gain with promising findings both in theory jU [5j [6] and in experiments [7J, El E]- 
Quantum direct detection [lQl [TTJ [12] by THz photon induced reduction in current 
has been experimentally observed [T3| [HI [15] and then applied to ultrafast detection 
and autocorrelation of short THz pulses[16]. A further topic on SSLs is strongly 
nonlinear dynamical phenomena that arise if one considers for example circuit models 
of SSL devices that incorporate effects that come from the interaction of the charge 
carriers. This gives rise to novel effects such as dissipative chaosfTTJ [18], and 
spontaneous generation of quantized [T9], EQl EH E21 [23], fractionally quantized, and 
non-quantized [221 121] DC bias. This spontaneous generation of a DC field is also the 
topic of this paper. 

Spontaneous generation of DC bias is the dynamical effect whereby a direct electric 
field appears due to a pure AC excitation. A significant point is that this spontaneous 
DC bias can in SSLs be nearly quantized, meaning that the total DC electric field, E& c 
follows the relationship 

fko , , 

E dc = T]n — , (1) 
ae 

where n is an integer, u is the frequency of the incident radiation field and the coefficient 
7] represents the fractional deviation from perfect quantization, rj = 1. The remaining 
constants are a, h, and e, which are the superlattice period, the reduced Planck constant, 
and the electron charge, respectively. This phenomena bears close resemblance to 
the inverse AC-Josephson effect [251 [26] . and is in fact just one example of similarity 
between Josephson junctions and SSLs[271[lEJEE]. Unquantized spontaneous DC has also 
been predicted for lateral semiconductor superlattices [291 EH] that only differ by their 
geometry from the bulk SSLs studied here. Similar effect has been earlier described in 
some models of homogeneous bulk semiconductors [3D] . 

We emphasize that this type of rectification occurs in perfectly symmetric 
structures, unlike classical rectification that relies on transport asymmetry due to 
contact charge inhomogeneities or nonlinearity. Rather, here a static DC field is 
formed dynamically due to presence of absolute negative conductivity (anc) at zero 
DC bias [HI 120] . anc in itself, is a well-known effect in SSLs [311 [32], [33], and has been 
theoretically found to be a robust effect in SSLs, occurring for a variety of miniband 
dispersions [M] and external excitations [351 1361 [37] . 

In SSLs the relationship (PQ) can persist for a wide range of parameters and so 
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Figure 1. Schematic figure of the SSL rectifier and the antenna. Incident radiation 
creates an oscillating potential C/ ex t and a fixed current idc is supplied via the leads 
connected to the antenna. 



this effect suggest the use of them as new type of frequency-to-voltage converters. A 
major hurdle exists, however, in the practical application of rectification in ssls: The 
parameters for which this dynamical rectification is expected to happen overlap exactly 
with the regions where the homogeneous spatial charge distribution is unstable. In this 
paper, we explore a possible way of avoiding this problem, and show that it should 
indeed be possible to observe the generation of a nearly quantized DC field with realistic 
SSL parameters. 

2. Effective circuit model with a dc current source 

As our model, we use an effective circuit calculation similar to [TU EHJ [22] where the 
active part, the SSL, is connected to an AC circuit modeling an antenna that couples 
the radiation field to the device. A schematic figure of the setup is presented in figure 
[TJ The external electric field generates an oscillating contribution {7 ext — U\ cos cut to 
the total potential across the device: U = U cxt — £E, where E is the electric field in 
the active region and I is its length. We take the amplitude V\ and frequency u of the 
external potential as our given parameters. Contributions from displacement currents 
are taken into account. Our new idea is to include a circuit that acts as a direct current 
source, providing a fixed current /dc- The DC voltage will then be determined by a self- 
consistent calculation of the model equations. We consider this voltage as the output 
of the device. 

The active part of the THz-rectifier is modeled using the well-known superlattice 
balance equations for sinusoidal miniband[3j 120] : 



The variables V and W are the average electron velocity and energy, respectively, where 
the averaging is done over a distribution function obeying the Boltzmann transport 
equation. W eq is the equilibrium energy, W eCL = -(A/2)I 1 (A/(2k h T))/I (A/(2k h T)), 
where A is the miniband width, T temperature, Boltzmann constant, and io,i are 




(2a) 




1 



(2b) 
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modified Bessel functions. E(t) is the total electric field that the miniband electrons 
experience, and it is obtained using the equation for current continuity across the active 
part of SSL: 

Above, Eo is the average DC dielectric constant of the superlattice, j is the conduction 
electron current density across the superlattice, R is a resistance that models ohmic 
losses in the superlattice, A is the area of superlattice cross section, and j ext is the 
external current density. This in turn is determined by considering the effective circuit 
model, which gives Aj ext = CJJ + I& c with C a being an effective antenna capacitance. 
The current density j is given by j = enV, where n is the electron density of the 
superlattice. Together, (jSJ) and ([3]) give a closed set of first order nonlinear ordinary 
differential equations. 

To facilitate analysis, we scale the variables V and W by their maximum values, 
giving W = (A/2)(l + w), V = a/S.{2h)~ l v. The new variables, scaled electron velocity 
and energy, v and w have the range — 1 . . . + 1. Electric fields are converted to frequency 
units by multiplying by ea/h: eaE/h = u. We obtain 

v = — uw — 7 v t>, (4a) 

w = uv - 7 w (to - w eq ), (4b) 

ii = — r 2 v — au + (f> ext , (4c) 

cxt = u s cos ut + i t, (4d) 

where r 2 = 0^/(1 + C a /C s ), the drive amplitude u s = eaU\/ {hl{l + C s /C a )), the 
electric field damping a = i? _1 (C s + C a ) _1 , the current bias io = aeId c /(M(C s + C a )), 
the superlattice capacitance C s = EoA/(4:7i£) and 

2ire 2 na 2 A 
h 2 e 

Here, u p \ is the plasma frequency and it describes the free oscillation frequency of 
the electron gas, and also serves as the parameter that determines the amount of 
nonlinear ity Note, that excluding the DC bias current term, equations (jlj) are formally 
equivalent to those used in [TTJ El] . For our simulations we use parameter values that 
closely match those of the superlattice studied in (39]: A = 120 meV, superlattice period 
a = 46 A, number of layers iV = 130, n = 9 ■ 10 17 cm -3 , 7 V = 15 THz, 7 W = 3 THz, 
and temperature T = 300 K. In addition we choose C a = C s ~ 5 fF, R = 1 kfl, roughly 
matching values used in [22]. Denoting r = 1/*/7v7wj the above values give rr ~ 2.25, 
ar ~ 0.015, and w eq ~ —0.75. 

In this paper we then consider the DC part u of the total electric field u, u = (u). 
We would like to find spontaneous generation of DC bias, i.e. (u) ^ 0, preferably nearly 
quantized so that (u) ~ nu. In the absence of the DC drive current, (H]) remains invariant 
in the transformation S: 

(t, v , w, u) — > (t + T/2, — v, w, — u), (6) 



u pI = — — • ( 5 ) 
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where T = 2ir/u is the period of the external AC field. If the governing equations remain 
invariant in transformation S, we say that they are symmetric or have symmetry S. In 
our previous works we have focused on the appearance of a spontaneous DC voltage via 
dynamical symmetry breaking^ZQ [TH [29j HDJ [28]. Let Uq be the generated DC voltage: 
u = (u) = T _1 f Q u(t') dt' (we use (■) to denote time averages over the period of the 
function being averaged). If the equations have symmetry S, a state with a nonzero 
average voltage, Uq ^ 0, implies that the solutions themselves do not follow symmetry S: 
{y{t + T/2),w{t + T/2),u{t + T/2)) ^ (-v(t),w(t),-u(t)). This is dynamical symmetry 
breaking. Here however, symmetry S is not present in the equations if i 7^ 0, and 
therefore one cannot speak of spontaneous DC following the breaking of symmetry. In 
spite of this, we will show that a key feature of spontaneous generation of a DC field 
remains: a sharp change in the total DC field near the conditions of symmetry breaking 
and also the quantization of that DC component. It should be noted that the DC 
bias current will only be essential during the turn-on of the device. Nothing prohibits 
disconnecting the bias current and restoring the symmetry of the problem after the 
device has settled to a steady state. 

The present work specifically addresses the problem of domain instabilities. An 
underlying assumption in the models used to predict the spontaneous DC field is that 
the electron distribution is spatially homogeneous. It is, however, well known that in the 
presence of negative differential conductivity (ndc), the homogeneity of the distribution 
tends to be violated and domains of different electric field strengths formfiTl W2\ H3] . 
This issue has been mainly ignored in previous works. For stable operation of the device 

THz-rectifier we then require that 

> 0, 7 

du dc 

where v^ c = (v) is the current averaged over its temporal period and u^ c is a DC probe 
field introduced to the total field u: u — Y u + u dc . We will refer to the steady state 
values of the realized DC electric field and current as the operating point. 

Condition ([7]) also effectively determines the appearance of a spontaneous DC field. 
Suppose v, w is a solution to the two first equations of (J4j) for a given field u containing 
a DC probe field represented by a slow function w dc . Substituting to (jlcl) one gets 
ttdc = — r 2 (v)(udc) — ctUdc- An equilibrium value of Ud c is stable if the derivative of the 
right-hand side with respect to Ud c is negative: 

- r 2 ^(u dc ) - a < 0. (8) 

OU dc 

Above, (v) is the time average of v, which is simply the DC current. Presence of the 
electric field damping a has the effect of stabilizing the DC field, however in the a — > 
limit equations (jHJ) and (j7]) coincide. Therefore, in order to achieve spontaneous DC 
from pure AC excitation, the superlattice needs to be driven into ndc. In the following 
sections, however, we will show that a small DC bias current can make the current- 
voltage curve (iv-curve) slope positive across a range of parameters while at the same 
time allowing for a substantial DC voltage to appear. 
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As the final point regarding stability against domains, we note that we should 
strengthen condition of (J7J). In an experimental setup the external AC field does not 
rise instantly, and so during turn-on and the initial transient phase the device might 
not be stable against domains. In the absence of a model that would take into account 
both formation of domains and spontaneous DC, we resort to comparing characteristic 
time scales for both processes to estimate which will dominate. Switching times r sw , 
that is, the time it takes for a spontaneous DC to appear should be faster than the 
domain formation times Td om to avoid the interference of the domains. In [20] where a 
similar model of a SSL rectifier was used, r sw was calculated to be at least in the range 
of hundreds of AC drive cycles T = 2tt/u in the case of heavily doped wide miniband 
superlattices. On the other hand, for similar superlattices the domain formation times 
can be as short as fractions of picoseconds [44]. Both of these characteristic times scale 
down as co p \ is increased, and since it is the dressed plasma frequency r, r < u p \, that 
here determines the switching time, we may suppose that r dom < r sw . This suggests that 
domains are likely to occur instead of spontaneous DC. This problem of ndc during 
turn-on and transient can be solved by assuming that the external potential is turned on 
slowly, and requiring that the device is stable for all drive amplitudes below the intended 
target value. This additional criterion bypasses all the uncertainty related to transients 
and should guarantee us stable operation. Next, we will show that in the limit u p \T < 1 
these requirements cannot be satisfied if a pure AC drive is considered, but can be if a 
small DC bias current is introduced. 

3. Weak DC biased switching in the uj p \t <C 1 limit 

First, let us consider the case where we only include the first harmonic of the field u. 
We assume the electric field has the form u = uq + u\ cos ut, and neglect the dependence 
of Ui on our main parameters u s and u. Also, we set 7 V = 7 W = 7- The first harmonic 
approximation is valid when the nonlinearity is weak, that is, when u p \t <C 1. 
Consider the DC part of the total electric field u. Averaging ( I4"c| yields 



where vq{uq,ui, lo) is the DC current under AC field u\ cosut. The resulting net voltage 
and the corresponding current are found at the intersection of the iv-curve and the line 



r 2 v oOo, Mi, w) = -au + io, 



(9) 



v r (u ) = r 2 (-au + i ). 
For a given electric field u, v (uq,Ui,u) can be explicitly written out as[3l 145] 



(10) 



OO 




(11) 



n=— oo 



v ET (u) 



= \w, 




(12) 
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The dependence vet(u) is the Esaki-Tsu current voltage curve [1] and in ( ITT]) we see 
the familiar form of an AC irradiated superlattice where the nth term represents an 
n-photon- assisted replica of vet[31 HEJ I4"T] . 

The effect of the DC bias is to shift the operating point to high DC field values as 
the AC amplitude is increased. This can be seen by considering a small DC voltage uq. 
Expanding ( II ip to leading order in Uq gives 

v (u , u u u) ~ goiux, u)u , (13) 

where go is the weak DC field conductivity of the SSL under DC and AC drives 

OO 

g (u 1 ,u) = J n (^"j v' ET (nu) (14) 

n=— oo 

and the prime symbol means the derivative with respect to u. In the weak DC field 
limit, stability condition ([7]) is equivalent to g (ui,u) > 0. Therefore, we wish to see 
how u changes as go tends to zero. The generated voltage m can be solved from (J9]): 

^0 = 2 7 l ° v , • (15) 

Setting go to zero one finds that u Q — > io/a, which is no longer a small quantity even for 
a relatively small i , since a C 1. This suggests that before the onset of ndc, the net 
DC electric field tends to large values. This is the basis of the DC biased slow switching 
scheme that we propose: The DC bias results in a large generated DC field that appears 
before NDC, thus avoiding the domain instability where go{u\, oS) < 0. Clearly, this alone 
is not enough, since the actual realized operating point may in fact sit on the ndc part 
of the iv-curve. Whether or not this is the case will be determined by ([7]) and (|9]). 

Let us now consider the frequency dependence of Uo- In the high frequency lot 3> 1 
limit the DC voltage can be very nearly quantized when the bias current is absent [20] 122] . 
This does not change when the DC drive is included. Since near uq = nu the n 
photon term dominates the series expression ffTT]) . the current is approximately given 
by J^(ui/u)vet{uo — nu). Using (]9]) and assuming |w — nu\ < 1 we find that 

io-anu 

M ~nw + 7 — — j-^f^- (16) 

cry + \w eq \J n (ui/ujyr z 

The correction term is small since we take ioT 2 ,ar 1, and so the realized operating 
point voltage is nearly quantized. Further, this point is always stable provided J n (u%/ cj) 
is not too close to zero. Note also, that ( 116]) gives an approximate optimal drive 
current for observing quantized Uo- The correction term vanishes and the quantization 
is improved by choosing i = i opt where 

z opt = nua. (17) 

The high-frequency limit is also particularly susceptible to domain instabilities during 
turn-on. This is because of hysteretic response to changing drive amplitude. The current 
replicas, that is, the individual terms in (|TT]) . define local peaks in the iv-curve. Since 
the slope of v T is low, the intersections between v T and vq tend to come in pairs, one on 
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Figure 2. [Colour online] Generated DC voltage as a function of AC field amplitude 
u\ and frequency uj. Remaining parameters are rr = 1, ar = 0.015, and w eq = —0.75. 
Regions shaded gray show stable operation. Dashed lines indicate borders of ndc 
[light red regions indicate NDC and domains], and solid thick lines show points of 
SN-bifurcations [dark red regions mark unstable branches of multistable operation]. 
Contours arc plotted at level values of u^/uj and are spaced at intervals of 0.1. (a) 
Bias current off, io = 0. Two regions of spontaneous DC are visible, however, they 
are flanked by regions unstable operation by all sides, (b) Bias current on, io = i opt . 
Strong response in uq is seen where spontaneous DC appeared in the unbiased case. 
Arrow at ui/ui ~ 1.70, ojt s» 1.62, uo/u f=s 0.39 marks the point where NDC first 
appears. Stable operation with large uq is then possible for lot < 1.6. 



each side of a local maximum of current. These intersections define a pair of operating 
points, one stable and the other unstable. Consider one such pair near uq = kcj. Should 
Jk{u\/u) tend to zero as U\ increases, that local current peak dips below v r and the 
operating points collide and are subsequently destroyed. This is essentially a saddle-node 
(sn) bifurcation and causes a problem since SN bifurcations always occur in conjunction 
with NDC. One can see this easily by noting that SN-bifurcation implies a double root 
of (EJ), which in turn means that slopes of Vq and v z must there agree: 
dv dc a 



du Ac , r 2 ' 



li 



'cic 

Clearly then, at an SN-bifurcation the system is already in NDC. Since equilibria 
are destroyed and created at SN-bifurcations, their presence implies multistability and 
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hysteretic response. Requiring stability against formation of domains therefore means 
that the voltage u response to changing u\ must be non-hysteretic. This effectively sets 
an upper limit to the frequencies by which stable turn-on can be guaranteed. 

In the low-frequency limit, lot <C 1, the current adiabatically follows the 
unirradiated iv-curve ( 1T21) . The DC current is then approximately given by (t>ET(wo + 
ui cos cut)). In this case, u as given by @ is always nearly zero, and exactly zero in the 
absence of the bias current. The low and high frequency limits together demonstrate 
that there is a band of frequencies near wr ~ 1 whereby strong response of uq can be 
observed without domains. In accordance with ({TBI) , the value of Uq should be closest 
to being quantized in the high end of that range of frequencies. 

To see what effect the DC bias current has on the parameters where rectification 
and quantized electric fields appear, we have in figure [2] plotted the DC voltage obtained 
by solving |9} as a function of u\ and u. We have chosen the parameter values as 
rr = l,u> cq = — 0.75,«r = 0.015. Borders of regions where ndc appears are indicated 
by dashed lines. Thick black curves mark SN-bifurcations. In figure EJja) the bias current 
is off. Two regions where spontaneous DC, \uq\ > 0, can be seen near u\/uj ~ 2.4 and 
ui/u ~ 5.5, the two first roots of Bessel J . The spontaneous DC forms plateaus 
satisfying the near-quantization condition (JTJ) in the high frequency end of the figure. 
The problem of domain instabilities is also evident. The regions of spontaneous DC are 
surrounded by regions of ndc. For low but non-zero values of u ndc is still present, 
but disappears for u closer to uj. Although the nearly quantized values of u are stable, 
they can only be reached via a hard-mode excitation, i.e. by instant turn-on and an 
initial electric field that is already close to the quantized value. With the more realistic 
adiabatic turn-on scheme, spontaneous DC appears impossible. 

Turning on the DC bias current essentially different behaviour is seen. This is shown 
in figure |2]^b) where DC bias current is set to i = i opt (optimizing for the u = u plateau, 
i.e. setting n — 1 in ( |T7|) ). The plateaus of nearly quantized spontaneous DC are still 
present where oj s /oj is close to the first two Bessel Jq roots. Additional plateaus closer 
to uo = 2u have now also appeared. The DC electric field is non-zero for other values of 
the parameters as well, but nonetheless there is a sharp increase or decrease in uq where 
spontaneous DC appeared or disappeared in the pure AC case above. Importantly, this 
time ndc is suppressed and stable operation is possible for a range of o; s from to 2.4, 
where the dependence u ~ u can be seen, ndc first appears at frequency wr « 1.62 
(indicated by an arrow in the figure, u\ ~ 1.70a;) followed by multistability at ur w 1.80, 
Ui m 1.68a;. This region extends long into higher frequencies and so it impedes reaching 
the quantized voltage plateau. Restricting to ojt < 1.6 the quantization is not so good, 
since well-developed plateaus are not present. Avoiding domains does then clearly limit 
how good quantization can be reached. 
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Figure 3. Generated DC voltage as a function of the drive amplitude at different bias 
currents. Solid curves indicate stable values while dashed curves operating points with 
ndc. Inset axis are the same as in main figure. Main figure: Parameters correspond to 
existing highly doped, wide miniband SSLs: lot = 2.0, rr = 2.25, j v t = 2.0, 7 w r = 0.5, 
ar = 0.015, and w cq = —0.75. A range of bias currents is from io = to io = O.I27 2 . 
Inset: Parameters are pushed to higher but reasonable values, lot = 3.0, rr = 5.0, 
7 v t = 2.0, 7 w r = 0.5, ar = 0.005, and w cq = —0.75. Bias currents range from io = 
to 0.397 2 . 

4. Numerical simulation, lo p \t > 1 

We now turn to the case of lo p \t > 1. Now all harmonics of the electric field are included 
and our results are parametrized by the applied AC potential amplitude o; s and frequency 
to. The problem is studied by numerically solving (j3J) for varying 10 and i§. We would 
like to find conditions for strong response in the DC voltage without NDC or hysteresis, 
and further, preferably quantized values of Uq. These requirements translate to having 
Mo dependence on u s resembling a step function. Results of the u p \t < 1 limit showed 
that quantization becomes better as frequency is increased, but that there is an upper 
limit in lot that still yields stable operation. We will use that finding as a guide and 
aim to find the highest applied field frequencies that are free from NDC and hysteresis. 
The cut upper limits turn out to be higher in the lj p \t > 1 case. 

We will start by considering a highly doped, wide miniband superlattice that is 
currently realizable. In figure [3] sample runs are presented where nearly quantized DC 
voltage can be achieved for parameters corresponding to such an SSL device. In the 
main figure, the values were chosen as rr = 2.25, lot = 2.0, 7 v r = 2.0, 7 w r = 0.5, 
ar = 0.015, and w eq = —0.75. A range of DC currents i from to O.I27 2 was used 
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to demonstrate the effect of current biasing. The io = curve shows again a region 
of spontaneous DC near u s /u ~ 2.4. This is surrounded as expected by ndc and here 
further by multistability and so stable operation is not available. As in the lo p \t < 1 
case, the inclusion of the DC bias removes the ndc and the multistability. For small 
but non-zero io the ndc is not completely gone, but does disappear quickly as the bias 
is increased. A rapid rise in w is observed near the symmetry breaking bifurcation of 
the ^o = curve (uj s /uj ~ 1.8), and further, this sudden increase quickly shoulders off 
to form a slightly tilted plateau where uq ~ oj. The realized uo depends only weakly on 
io, showing that the device is not sensitive to the choice of bias current as long as it is 
sufficiently high to remove the instabilities. Higher near-quantized plateaus for which 
Uq/u ~ n with n > 1 also become available, though now the quantization is less good. 

The inset in figure |3] shows what can be achieved with a reasonable push in the 
parameters. Here, we have boosted the plasma frequency to rr = 5 and drive frequency 
to lot = 3.0 while reducing the inverse parallel resistance a to 0.0057. Dependence on 
the bias current is again similar (here, io = . . . O.397 2 ) with the case of zero external 
current showing instabilities that are removed as io is increased. With these parameters 
the shoulder where uq levels off to ~ w is sharper and the quantization better, yet uo is 
single valued in co s as desired. Quantization with r\ 0.95 (u = r/u) can be observed 
here for a wide range of u s . 

Note that the u p \r < 1 limit suggests a scaling law. The equations that determine 
Uo, Q and ffTT|) . depend on r, a, and i only via the fractions a/r 2 and io/r 2 . Therefore, 
one could hypothesize that the increased stability of oj p \t > 1 case is only due to a 
decrease in the effective inverse parallel resistance a/r 2 , and could then be replicated 
by keeping r fixed and decreasing a instead. We tested this hypothesis by running 
simulations with fixed a/r 2 and varying r and found it to be false, so that increasing r 
really does have an effect of suppressing ndc. 

We see that stable operation at higher applied field frequencies can be obtained in 
the limit of a high plasma frequency, or a high doping and a wide miniband. Because 
higher lot can be reached, also the quantization of the spontaneous DC voltage is 
improved. Nonetheless, ndc and hysteretic response still set strict limits to frequency 
and hence to the quantization that can be achieved. Stability against domains must be 
balanced against how small \uo — nu\ is desired. To elaborate further on the parameters 
that affect the spontaneous DC, we note that there is some sensitivity to the value a. 
Zero inverse parallel resistance, a = 0, results in a more hysteretic response, and so 
a small degree of damping in the equation for u is preferable. However, for ar > 0.1 
the generated DC voltage becomes close to zero, so smallness of a is still required. In 
terms of the DC bias current, the device does not appear to be very sensitive to its value 
as long as it is sufficient to remove ndc near zero DC voltages. Increasing i Q further 
to moderate values has little effect on stability, it only shifts the point where strong 
response occurs to lower values of lo s /oj. Again, we note that since the spontaneous DC 
states are often stable when io = 0, the bias current is not necessary anymore once a 
stable equilibrium is reached provided of course that the unbiased operating point is not 
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in NDC. 

5. Conclusions 

We have presented a turn-on scheme for a semiconductor superlattice THz-rectifier that 
addresses the problem of stability against formation of electric domains due to NDC. 
While a spontaneous DC voltage state of a pure AC driven semiconductor superlattice is 
stable, this voltage state can only be reached by driving the device into conditions where 
domains do form. Our proposed solution to this problem is to apply a weak DC current 
bias while the applied field is adiabatically increased to its target value. The external DC 
bias has the effect of shifting the operating point towards higher DC voltages whereby the 
device shows positive differential conductivity, and therefore is stable against formation 
of domains. This shifting appears as a strong response in the DC voltage as a function 
of the amplitude of the external AC field and is seen for parameters that correspond to 
dynamical symmetry breaking in the pure AC driven case. 

For typical superlattices and experimentally realizable parameters of applied field, 
together with small DC bias currents, the strong response in the DC voltage and stable 
switching to a spontaneous DC field state are indeed found. Appearance of NDC and 
hysteretic response however do still set an upper limit to the applied radiation (angular) 
frequencies uj that give the desired behaviour. For weak nonlinearity, u is roughly limited 
by 1.5 ~ 2 times the characteristic scattering rates. Somewhat counter-intuitively, 
in the limit of strong nonlinearity (high plasma frequencies) stable operation can be 
achieved for higher drive frequencies. In accordance with previous results on generation 
of spontaneous DC voltage, we found the generated DC field becomes closer to an integer 
multiple of hu/ea as u> is increased. Therefore, our results show that for optimal, 
stable frequency to voltage conversion, high plasma frequencies are preferable. Our 
numerical results show that good quantization of the generated voltage can be achieved 
for realistic semiconductor superlattices. These findings demonstrate that it may indeed 
be possible to apply semiconductor superlattices as room temperature frequency-to- 
voltage converters and detectors of a strong THz radiation. 
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